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^-H . Abstract 

o, 

^ ! We study the stability of the ordinary Landau Fermi liquid 

phase for interacting, spinless electrons. We require causality 

rH ; and demand that the Pauli principle be obeyed. We find a 

JL ; phase diagram determined by two parameters: the particle 

Ch ! density and the interaction strength. We find that the ho- 

CJ ■ mogeneous, constant density Fermi liquid phase of a spinless 

Fermion gas is never stable in ID, but that it may have a 

rS ' restricted domain of stability in 2D. 
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INTRODUCTION The theory of the homogeneous electron gas, the "Fermi liq- 
uid" , plays a very important role in study of metallic properties in real solids. For 
even within the continuous jellium model, Wigner [1] argued long ago that at low den- 
sity, the translational symmetry of this phase is broken and electrons are localized in 
a lattice to lower the potential energy. Qualitatively speaking, as the mean interpar- 
ticle distance r^ oc 1/n^/'^ increases, i.e. density n decreases, the Coulomb potential 
energy (oc 1/ts) of such an electron gas will eventually dominate the kinetic energy 
(oc 1/rf). As a result, the homogeneous Fermi liquid phase with a well-defined Fermi 
surface is unstable in the low density limit. The nature of the phase which replaces 
it is, however, a matter of conjecture. The localized particles in a Wigner lattice 
are one possibility. A two-phase regime (high density/low density) is another. The 
marginal Fermi liquid [2] and the quantum Fermi liquid [3] in addition to a variety 
of "anyons" are distinct possibilities in low dimensions. In 3D, there have been a 
number of estimates of the critical r^, as summarized by Care and March [4]. 

Sometime ago Mattis [5] proposed a lattice model which incorporates the "gran- 
ular" nature of electronic states, reflecting the finite size of the Wannier orbitals of 
the conduction band of a simple metal. He found that Coulomb interactions alone 
could lead to instability of ordinary paramagnetic Fermi liquid against various inter- 
esting phases, including charge-density waves. Following a theorem of Ginzburg and 
Kirzhnits [6], he proposed a loose but general stability criterion, viz.: the effective 
two-body interaction v{q) is required to be non-negative for all wave vector q. 

In this paper, we propose a refinement taking into account the short-range corre- 
lation in addition to the causality requirement considered previously. The application 
is to spinless electrons living on a simple rf-dimensional cubic structure undergoing 



mutual Coulomb interactions. In contrast to the jellium model where the only length 
scale is r^, there appears an additional parameter in the lattice length a. (In the case 
of real electrons with spins, the magnetic degrees of freedom interfere with the charge 
dynamics to make the analysis much more complicated.) 

The stability criterion of a spinless Fermi liquid will be based on simultaneously 
satisfying two requirements: (a) by causality, the inverse static dielectric function 
l/£:(g, 0) < 1; (b) by the Pauli principle, the zero-range pair correlation function of 
spinless Fermions 5'(0) = 0. Condition (a), the "Ginzburg-Kirzhnits criterion" [6], is 
proved using the Kramers- Kronig relation: 

where "P." stands for principal value. As Im{e(g,a;)} > 0, it follows that the recip- 
rocal static dielectric function l/e(g, 0) < 1. 

The function g{r) is defined [7] as the probability of finding two particles separated 
by a distance r, normalized to be 1 asymptotically at large r. Fortunately g{r) can 
also be determined by the dielectric function according to a well-known many-body 
identity (see Eq.(1.6.2), (5.4.12) in reference [7]): 

9{r) = l + U^Y.^''^-'S{q)-l) (2) 

where 

^, , f^ du 1^,-1-, /^^ 
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and v{q) is the two-body vertex function. 

It is well-known that condition (b) is never satisfied within the jellium model 
at low density. The calculated [7,8] pair correlation function g{r) of the continuum 



theory inevitably turns out to be negative at small r, evidence of an unphysical 
excess of exchange and correlation holes around each electron. However, this defect 
is easily cured when electrons live on a discrete lattice where the short-range, on-site 
interaction and the long-range, off-site Coulomb interaction can be treated separately. 
For we have observed that within RPA, condition (b) is equivalent to the requirement 
that the ground state energy be stationary with respect to a fictitious zero-range 
interaction parametrized by a new coupling constant a. This appears to cure a major 
flaw in the RPA, allowing it to be used for semi-quantitative purposes. 
SPINLESS FERMIONS The extended Hamiltonian is: 

^ = Y^{^k- ^)clck + 2 H UijUiUj + -J2ani{ni-1) (4) 

k iy^j i 

The third term X]i^i(^j ~ 1) = vanishes for spinless Fermions and a is the cor- 
responding undetermined Lagrange multiplier. The c's are Fermion operators and 
rii = c\ci is the particle number operator on site i. In ID, e^ = —2t cos ka (a the lat- 
tice length) and t is the transfer integral. The chemical potential fi = ep, the Fermi 
level. The two-body Coulomb interaction Uij = e^/vij, rij = \i — j\a; generalization 
to higher dimensional hypercubic lattices is straightforward. 

To find the ground state energy, we rewrite H in momentum space keeping careful 
track of all the terms: 

H = J2ie, - ^ylc + ^ ^(a + A\/(g))p(g)p(-g) - ]-aNn{l - n) (5) 

k qj^O 

where the generalized Ewald sum 

A\/(g) = ^e*«-"'— (6) 

and charge density operator p{q) = Z^fccl+^Cfc. A^e is the total particle number, A^ 
the total number of cells and n = N^/N is the particle density. In the usual way, we 



estimate the correlation function 

1- [a + AV[q))Il[q,u) 
within the RPA approximation by taking 11 ^ 11'^, the Lindhard polarization function 
[7]. Then the ground state energy is given by an integral over the coupling constant 
yielding: 

r°° duo 

with Eq the energy of free Fermions. Optimizing Eq with respect to a, 
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EG = Eo + -lmJ2 — log(l - (« + AV{q))U\q,uj)) - -aNn{l - n) (8) 
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we find: 



lj2Sa{q) = l-n (10) 



where 



Q ( )= H— 1 T zl nn 

"^^^ Jo TT ■n(a + A\/(g)) "^l-(a + A\/(g))no(g,cu) ^ ' 
is the structure function. The dielectric function which appears in Eqs.(l)-(3) is: 

e„(g,cu) = l-(« + A\/(g))n°(g,cu), in RPA. (12) 

By combining Eqs.(2) and (10) we find: 

^«(o) = i + ^(^E%)-i) = o (13) 

^ ^^ q^O 

Therefore 

_^ = ^ ^a(0) = (14) 



and the stationary condition enforces the Pauh principle. Using the value of a which 
satisfies this criterion, we now identify Vaifl) = a + A\/(g) as the "best" effective 
two-body vertex function. Since n'^(g, 0) < 0, Ginzburg and Kirzhnits' criterion 
5^(g, 0) = 1 - ■I4(g)n°(g, 0) = 1 + ■I4(g)|n°(g, 0)| > 1 implies \4(g) > for all q for 
the system in the Fermi liquid phase, just as in reference [5]. The interpretation of 
our result is that the "best" choice of a compensates for the errors in RPA at short 
distances, and transforms it into a more reliable procedure. 

NUMERICAL RESULTS The problem has thus been reduced to to find a solution 
for a which satisfies both 

K.(g) = a + ^V[(i) > 0, Vg, and ^^(0) = (15) 

We will first solve for a, then plug it into the inequality to check whether it is satisfied. 
If not, some symmetry breaking must occur. 

We present the numerical results in ID and 2D. In ID, the generalized Ewald 
summation is: 

„2 



^iq-na 



A\/(g) = y: 

~ll \na\ 

= m 4sm — 16 

a ' 2 ' ^ ^ 

and it has a minimum AV{'n-/a) = — 1.386e^/a. In addition, the retarded polarization 
function 11° (g, u) can also be evaluated analytically, using the tight-binding spectrum 
ek = — cos ka (taking 2t = 1 as our unit from now on). By symmetries n°(g, u;) = 
n*'(— g, cj), Ren°(g, — cj) = Ren°(g, cj), Imn°(g, — c<j) = —Imll^{q,uj), and because of 
periodicity in q, we only need to consider a; > and < qa < tt. For uj > 2 sin y'- 

Relf{q,uj) 
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(arctanf , ^ — arctanf , ^ ) (17) 

(cj/2sinf)2-l \/(cc;/2sinf)2-l 



Fortu = 2sinf : 
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27rsinf' cos(f + A;^) cos(f 



For <cj <2sinf : 



27rsinf Jl-(cu/2sinf)2 



Ren°(g,cu) 

cos(f - kp) - Jl - (cu/2 sin f )2 cos(f + A;^) + Jl " (^/2 sin f )2 

In I . , Kl9) 

cos(f - kp) + y^l-(cu/2sinf)2 cos(f + kp) + y^l - (cu/2sinf )2 

The imaginary part is: 
Imn°(g, a; 
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Of qa . qa i \q( i Qa . qa 

•y(cospcos h smpsm coskFJyycoskp — cospcos h smpsm — ) 

Zi Zj Zj Zi 

■^(2sin^-cu)| , / , , ,, . ,„,2 (20) 

In 2D, the Ewald sum AV(g^,gj^) has been given by Glasser [9]. But it was 
expedient to input a shghtly "screened" form of AV(g^, qy). This does not affect the 
physics but simplifies the numerical work: 

{±M,±M) 2 ^ 



(m,n)7^(0,0) 



i^m^ + n^ 



having a minimum Ay(7r,7r) = — 1.024e^/a for M = 10, q^ = 0.69 (the numbers 
adopted in our numerical simulations.) Similarly, the lattice Il^{qx, qy, uj) (a general- 
ized Watson integral) is replaced by the isotropic 2D Lindhard polarization function 



H^^Jq"^ + qy,ij) [10] evaluated in the effective mass approximation, and the filling 
factor is determined by n = kp/Arr for n < | (electrons) and 1 — kp/Air for n > | 



(holes.) These approximations should be good for all almost empty or almost full 
band, but should not be considered accurate at or near half-filling. 

Only two parameters govern our problem: one is the particle density n. The other 
is U /W = e^/2ta, measuring the Coulomb interaction strength U relative to half 
band width W. 

The solution a is presented in Fig.l for ID case, and Fig.2 for 2D case for n < ■^. 
The results are symmetric about n = |, so n > | is not shown. 

As we see in Fig.l, solutions a's never reach 1.386, the maximum of —AV in ID. 
This implies that causality can not be satisfied as there always exist some Qc < vr/a 
such that for q > Qc, V{q) < for all density n and interaction U/W. Thus we 
see that the ID spinless homogeneous Fermion gas is unstable. This instability is 
already known by way of the exact solution by Mattis and Lieb [11] of the Luttinger 
model. At arbitrarily small values of the interaction energy, these authors found that 
the discontinuity of < n^ > at the Fermi surface disappears, and the Fermi liquid 
is transmuted into a "Luttinger liquid" with a corresponding charge in correlation 
functions. 

Fig.2 shows the situation in 2D; the a's increase along with n, and finally exceed 
1.024, the maximum of —AV in our calculations, after which they satisfy causality. 
Hence, at a given interaction strength U/W, the spinless Fermion gas is unstable 
only at density n lower than some critical ric- The phase diagram in Fig. 3 shows 
Uc v.s. U/W. As U/W increases, Uc approaches |. An additional possible region of 
instability against dimerization near n = | can not be investigated within the effective 
mass approximation. Thus it is possible that at sufficiently large U/W, the 2D Fermi 
liquid becomes altogether unstable, but due to the errors in obtaining the optimum 



a at large n using the effective mass approximation, and because of the systematic 
and inherent errors in the RPA at large U jW ^ we can only speculate on the strong- 
coupled behavior. Efros [12], who has studied the strong-coupling limit of the 2D 
electron liquid (in a strong magnetic field which effectively quenches the spin degree 
of freedom) has found that in the absence of motional energy, the electron-electron 
interaction leads to a very inhomogeneous phase. His result may be complementary 
to our RPA approach, which was predicated on the weak-coupling regime. 

The effects of introducing a local field correction [8] G'(g) were also investigated, 
but were found to have no qualitative infiuence on our results. 

In conclusion, we find for spinless Fermions that the Landau Fermi-Liquid phase, 
characterized by a one-to-one correspondence with the Fermi gas, including the dis- 
continuity in < n^ > at kp^ has only a restricted region of validity in 2D, and that 
it is never applicable in ID. This last result was already anticipated in the solution 
of the ID Luttinger liquid [11], albeit from a different point of view. In the present 
paper, we reach these conclusions by optimizing RPA using a Lagrange multiplier a. 
At the optimum a, the RPA is found to satisfy the Pauli principle, thereby becoming 
a simple and satisfactory tool with which to study the stability of the Fermi liquid. 

The phase diagrams for SU(2) Fermions in 2D and in 3D are more complex and 
require a subtler analysis. They will be examined separately. 
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Figure Captions 

FIG.l 
Solution a (in unit of e^/a) required to satisfy QaiS^) = in ID. Solid line is for 
U /W = 1.0, dash line for U/W = 0.1. Causality requires a > 1.386, which can never 
be satisfied. 

FIG.2 
Solution a (in unit of e^/a) for (?a(0) = in 2D. Solid line is for U/W = 1.0, dash 
line for U/W = 0.1. Causality requires a > 1.024, as shown by horizontal dash line. 
It fails to be satisfied at low density. 

FIG.3 
Phase diagram in 2D, indicating maximum region of stability of Fermi liquid (there 
may be an additional dimerization instability near n = ^ which we are unable to 
examine by our methods.) 
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